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ABSTRACT 

Conventional data obtained in 1932 are contrasted with SEASAT-SASS and SMMR 
data to show how observations at a single station can be extended to an area of 
about 150,000 square km by means of remotely sensed data obtained in nine 
minutes. Superobservations at a one degree resolution for the vector winds are 
estimated along with estimates of their standard deviations. From these super- 
observations, the horizontal divergence, vector wind stress, and the curl of the 
wind stress can be found. 


INTRODUCTION 

The basic tenets of those who believe that it is possible to predict the 
weather are three. These are (1) that the time evolution of *i\e meteorological 
variables can be correctly described by appropriate equations, (2) that the 
initial conditions for the start of the prediction can be adequately specified 
and (3) that computers can be made with a large enough capacity and a fast 
enough speed to compute the evolution of the weather for the future. These 
tenets evolved with time from the days of Bjerknes (1904) who wrote that what 
was needed was "a sufficiently accurate knowledge of the state of the atmo- 
sphere at the initial time" and "a sufficiently accurate knowledge of the laws 
according to which one state of the atmosphere develops from another", through 
Richardson's (1922) "Weather Prediction by Numerical Process", on to the day 
when the authors of a recent paper wrote "The main reason we cannot tell the 
weather tomorrow is that we do not know the weather today". 

These tenets are interdependent because improved computers provide higher 
resolution. Improved descriptions of the physics of the atmosphere require 
more data. The problem is to predict the synoptic scale meteorological evolution 
and to filter out the higher wavenumbers and frequencies. This is a corrolary 
of the basic tenets and involves all aspects of turbulence, where much is yet to 
be learned. 

Much of what has been learned in the past about the behaviour of weather 
systems has been lost because of present day computer constraints. To our 
knowledge, no operational hemispheric numerical model with a one degree by one 
degree resolution, and a corresponding vertical resolution, presently exists. 
With an adequate description of the physics of the atmosphere and of subgrid 
turbulence, Including the possibility of transfer of energy from high to low 
wave numbers, a one degree, at least hemispheric model, is presently conceiv- 
able. 


OCTOBER 21, AND 22, 1932 AT HOLYHEAD, ENGLAND 

Before showing what can be achieved by means of remote sensing, it is useful 
to go back to a text by Brunt (1942). It is then possible to compare the results 
of the analysis of SEASAT data obtained in about 9 minutes with data obtained 
over a 16 hour period in England. 
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Fig. 1(a) Synoptic analysis for the weather on October 22, 1932 at 0700 for the 
British Isles and part of Europe and (b) time histories at Holyhead. Wind speed 
from 0 to 50 miles per hour; Direction, meteorological convention; Temperature, 
degrees Fahrenheit 50° to 58°; Integrated rainfall in millimeters, 0 to 10 mm; 

Sea level pressure, millibars, 990 to 1006. From Brunt (1942). 





Fig. 1(a) shows the synoptic scale analysis for 07h Oct. 22, 1932, London 
time, for an area around the British Isles. Fig. 1(b) shows the Lime history of 
the evolution of winds, temperature, integrated rainfall and sea surface atmo- 
spheric pressure at Holyhead, England (approximately 52.5°, 4.8° W) beginning on 
Oct. 21 at 18h and ending on Oct 22 at lOh (16 hours). Some data have probably 
been omitted for clarity in Fig. 1(a) but there are roughly four observations 
per 10 degree latitude by longitude "square” which permit an adequate analysis. 
The cold front of the open wave cyclone is In the process of passing Holyhead, 

The anemometer had a rather short time constant so that mesoscale fluctuations 
of the wind were recorded. Even in 1932, it would have been obvious that a two 
minute average could give erroneous synoptic scale values as recently rediscov- 
ered by Pierson, et al. (1980) and made quantitative by Pierson (1983a). 

If one imagines the wave cyclone to be moving without too much change in form 
along an appropriate track, the winds in advance of the warm sector were ENE for 
the first eight hours building from 6 or 7 miles per hour (2.24 mph equals lm/s) 
to 22 mph (with gusts to 29) followed by a decrease to 6 mph and a wind shift to 
the ESE. The wind shift to SSW occurred in about 20 minutes. The temperature 
rose slowly from about 1:40 to 3 AM then rapidly for 10 to 15 minutes with a 
1 - exp(-Bt) asymptotic approach to 56° F. The total warming was 6° F. 

The winds in the warm sector were very gusty and turbulent with 40° fluctua- 

tions in direction and speeds that varied moment by moment from under 10 mph to 
nearly 45 mph. The average wind in the warm sector increased from about 20 mph 
to 25 mph before the cold front passage. 

The cold front went by in two steps at about 6:30 and 7:20 with a wind shift 
to the WSW followed by the final one to the NW. The temperature dropped sharply 
twice in 2° F steps and settled down to a value of 52° F. 

Light pre-warm front rain occurred from 1900 to 1930 and from 2050 to 2255. 

Heavier rain began after midnight and lasted steadily for more than two hours 
tapering off after the warm front passage and stopping at about 0620 at the first 
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passage of the cold front. 

With modern communication systems, it would he possible, ip principle, to 
obtain similar data with data buoys over the ocean and closely spaced land 
installations for pre processing locally to define the planetary boundary layer 
with sufficient accuracy at the synoptic scale and to collect and iu>e r lu data 
in time for forecasts. The cost would be substantial. Over the oceans, remote 
sensing can provide a cost effective alternative. 


i 


! 
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THEOPY 

To predict the weather, meteorologists start with the hydrodynamic equations 
for a rotating Earth plus l^O conservation equations, sources and sinks, plus 
radiation budget equations, plus whatever else might be needed and try to put 
them in a form that can be integrated forward by means of about 10 minute time 
steps so as to reach out to the future. The first difficulty is that the hydro- 
dynamic equations are inappropriate. It is impossible to measure the various 
fields at the resolutions they need, and equally impossible to integrate them 
forward in time to make a forecast. The atmosphere is turbulent at scales that 
are, and always will be, unresolvable for such purposes. The equations that 
are actually used require time and space averaged values and therefore become 
strangely altered because they must account for the effects of fluctuations 
within the space and time average. A host of new terms such as -^u'w'^* arise 
because of nonlinearity that must be removed by means of an appropriate closure 
model. Just how this ought to be done is a matter for discussion as in Robinson 
(1978) , among others. 

To put things simply, the synoptic scale winds are needed for a numerical 
weather prediction as defined in a circular way by Pierson (1981a). Each SASS 
wind in a swath is defined by four numbers; they are the latitude, 0 ^ , longitude 
Ai, wind speed at 19.5 meters for an effective neutral wind |v|^, and the mete- 
orological wind direction, Xi* after the removal of ambiguities. 

Each SASS wind in a data set departs from the desired synoptic scale value at 
0^, Xi f° r two reasons (at least). The footprint (or cell) illuminated by the 
radar is analogous to a rather brief time average that would not quite recover 
the synoptic scale value. For a typical cell size and a synoptic scale wind of 
15 m/s the variability due to this effect can, for example, be ±0,38 m/s at one 
standard deviation. Other sources of error caused by the radar itself, and 
lumped into one term (Pierson(1983(a) , 1983(b)), produce about an 0.6 m/s varia- 
tion, but this can vary for various reasons as a function of the wind direction 
and incidence angle. For direction, the corresponding values are 1.2° and 6° 
for this example. The combined effect is about 0.71 m/s for speed and slightly 
over 6° for direction. These are representative conditions, but over a full 
swath these two effects can combine to give variations more than twice as large 
or less than half as much. They are inseparable after the vector wind is found. 
Neither effect is desirable for the purpose of synoptic scale analyses. Both 
are random from one measurement to another. 

The data analysis procedure reduces these sources of variability and recovers 
a reliable estimate (in a statistical sense) of the synoptic scale wind at 19.5 
m above the sea surface for a neutrally stratified atmosphere. For each integer 
value of latitude and longitude in the swath, say 0 O , > Q , all SASS wind within 
the 2 by 2 degree M square M surrounding the target point were found. The average 
latitude and longitude of the SASS winds in this sample would not quite be at 
6 0f X 0 and a few values were dropped from the sample at the outer edges to 
obtain an average latitude and longitude close to 0 O , X 0 . This slightly reduced 
sample was then processed as a superobservation to obtain synoptic scale fields 
with estimates of sampling variability effects at a one degree resolution by 
overlapping the 2 by 2 degree squares. 

For each sample a mean direction can be found as in (1). 
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( 1 ) 

to X such 

( 2 ) 

(3) 

and their variances can be found as in (4) and (5). Also quantities such as 
£(6^ - Qj = ^Ae^ , + TFT, A Q + AA, VAR (A0)> and so on, are needed because the 

(4) 

(5) 

From standard statistical theory, the expected value of the mean is the mean 
and the variance of the mean is the variance of the sample as in (4) and (5) 
times N”l. The only difficulty is the presence of gradients in the winds over 
the swath which can overestimate the variances that are needed and the slight 
mis-location of the desired winds by 7RF, ST. 

Further data processing in terms of east-west components, and north-south 
components, Vq , allows the effects of gradients to be removed to first order, 
the wiM vector to be located at X Q , 0 Q , and the standard deviations of the mean 
east-w<st and north-south components to be found. The final results, with added 
complexities at the edges of the available data, were values for the vector 
components of the wind at each 0 O , X Q that were much more accurate estimates of 
the synoptic scale wind with known effects of sampling variability as in equa- 
tions (6) and (7). (We have simplified notation for this paper). 
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In (6) and (7), U. and V A are the desired, but always unknown, synoptic 

>rv, es - 


scale values, U anc 


are their superobservation estimates and AU- , AU. 


AV^, express the reduced effects of the components of (4) and (5) after removing 
the contribution from gradients in the wind field and obtaining the appropriate 
vector components. The values of t^ and t£ near the middle of the swath where 
N can be above 30 are close to being values drawn at random from a zero mean 
unit variance normal probability density function. 

As an example, at 54° N 209° U (for another SEASAT pass), the superobservation 
yielded a wind from 260° at 20.2 m/s. The other data obtained then gave (8) and 
(9) in m/s. From (8) and (9) , as one of many examples that could be given, the 
uncertainty in wind speed is ±0.17 m/s at one standard deviation and the uncer- 
tainty in direction is ±2.75°. It will be a while before meteorologists will 
believe such results. 


U Xs - 19.95 + t 1 (0.17)+ t 2 (0.14) - 19.95 + t*(0.22) (8) 

V 0s - 0.156 + t 1 (0.02)+ t 2 (-0.96)«0.156 + t*(0.96) (9) 

By these methods as an example for the REV to be studied, the original 1,081 
SASS winds were reduced to 120 superobservation winds as in (6) and (7) with 
each wind accompanied by estimates of its sampling variability for a reduction 
to 11Z of the original data for further processing. 
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For every point in the super observation wind field that had four points 
around it, it was possible to estimate the horizontal field of divergence at 
19.5 meters and to carry along the four terms in (6) and (7) so as to obtain 
the standard deviation of the estimate of the divergence. The appropriate 
equation is (10) 
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The finite difference form is (11). 


div,( W (ex))- 
l n o o s 

2.5-10" 6 

'<V 9 o- 

cos 0 

■+• cos 0 (V (0 + 1 

•V S0 o 

, X ) -V (0 -1,X 

O S9 O. O 

- D , 0349065 Bin 0 

0 

W ; 

V 


1.(0 , X -1)) 
sX o* o 


(ID 

Five different terms such as (8) and (9) are needed. The expected value of 
the 10 differents t's that enter is zero so that the expected value of the 
divergence can be found. 

Also 
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can be found as a weighted sum of the squares of the AUj and AVj so that the 
standard deviation of the estimate of the divergence can also be estimated. 

The field of horizontal divergence is rarely calculated from conventional data 
because of large errors in measurement. The divergences that were found formed 
logical coherent patterns. 


THE WESTERN NORTH PACIFIC FOR SEPTEMBER 19, 1978 AT ABOUT 1800 GMT 

On this date, SEASAT crossed the Equator, on its 1212th revolution of the 
Earth in a retrogade circular orbit with a period of about 100 minutes at an 
altitude of 800 km., going north bound, a few seconds before 1754 GMT and ob- 
tained the data that were analysed for the following material between about 1803 
and 1812 GMT from the SASS and the SMMR. For details, see for example IEEE J. 
Oceania Eng. Vol. OE-5 and J. Geophye. toe. Vol. 87, No. C5 and Vol. 88, No. C3. 
Further details on these particular results and for three additional weather 
patterns are given by Pierson, et al. (1983). During this 9 minute interval, 
the SASS found 1,081 values for the winds over an area about 600 km wide ex- 
tending from 35° N to 60° N in a one side scanning mode. The SMMR also obtained 
data over a slightly narrower swath. At nearly the same time, the GOES space- 
cr *ft obtained both visible and infrared lmagry. The results for the synoptic 
scale wind field and the field of the horizontal divergence of the wind will be 
described. Fields for the wind stress, (u^)**, and the curl of the wind stress 
will be found by similar methods, also. 

Me Mur die and Katsaros (1983) have used the SMMR data to obtain the values in 
Fig. ^ for this REV. Their results can be combined with the vector wind values 
snd the divergence field to produce data over an area of roughly 150,000 square 
kilometers that are effectively the equivalent of having records similar to Fig. 
1(b) from every one degree latitude longitude Intersection in the swath for 
synoptlcally averaged winds and for tha rainfall ratas. Given the few ship 
reports of sea surfaca atmospheric pressure within the swath. It ought not to be 
unduly difficult to recover the sea surface pressure field within tha swath. 

Air temperature just above the sea surfaca is still a problem for remote sensing, 
and tha sparclty of conventional ship reports makes such data difficult to ob- 
tain. 
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Fig. 2. (a) Geostationary cloud image of the 

extratropical cyclone for SEASAT REV 1212. (b) 

Streamline, isotach analysis of the superobserva- 
tions. Speeds in meters per second. Data from the 
SEASAT A Scatterometer , (SASS). (c) Field of 
horizontal divergence computed from the superobser- 
vations. Fig. 2(a) is from Woiceshyn, et al. (1979). 


The GOES cloud imagry, the SASS superobserva- 
tion winds, and the divergence field computed 
therefrom plus the SMMR lesults provide a detailed 
meteorological description of an occluded cyclone 
approaching the vest coast of North America. The 
GOES cloud image in Fig. 2(a) shows an occluded 
low crossing the vest coast of North America. A 
long cold front ’’.rails to the south and west over 
the Pacific. The cloud shield extends well into 
Canada and covers Queen Charlotte Island and most 
of Vancouver. 

Streamlines and isotachs (m/s) are shown in 
Fig. 2(b). Strong warm air advectlon for the 
streamlines surrounded by the IS m/s isotach and 
to the south to 10 m/s is shown. The wind shift 
in the general vicinity of the occlusion is shown. 
* Queen Charlotte Island is not shown for simplicity. 

' The field of divergence is shown in Fig. 2(c) based on 65 values ss a composite 
with the frontal analysis based on conventional analysis plus a schematic of the 
cloud image shown in Fig. 2(a). Strong convergence is shown at the intersection 
_ of the cold and warm front with a band less than -4 X 10 -5 sec' 1 extending over 
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The estimate was between 1 X 10~5 
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Fig. 3. The integrated water vapor in k^/rn 
(a), the integrated lio’iid water in kg/in'- (b) , 
and the rainfall rate in mm/hr (c), found from 
the Scanning Multichannel Microwave Radiometer 
on SEASAT. In (a) the contours that parallel 
the cold front have values of 16, 20, 24 and 28. 
The dark area has values over 32 and the inner 
white area values over 40. In (b) , the area in 
the center has values of 0.8, 1.2 and 1.6 with 
1.6 in white. In (c) , the rainfall rate contours 
are for 0.4, 0.8, 1.2 and 1.6 nm/hr. with the 
center white dot showing 1.6 mm/hr. From 
McMurdie and Katsaxos (1983). 

more than 5° of longitude. Convergence is 
also shown to the west of Vancouver in the vara 
air sector of the system. Divergence values 
as high as (3.55 i 0.38) X 10* 5 sec~l and as 
low as (-7.98 ± 0.91) X 10* 5 sec-1 were obtained. 
The standard deviations show that tha firs'- was 
surely positive and the second surely negative. 
Values of near zero divergence would be indi- 
cated if the estimate plus one standard devia- 
tion and the estimate minus one standard devi- 
ation changed sign. This happened sir times. 
sec-1 and 1 X 10"5 8 ec-l 28 times, suggesting 


areas of nearly nondivergent flow. 

The higher convergence values would produce upward vertical velocities of 1.6 
cm/sec at a height of 200 m. Air moving at 15 m/s will rise 57 m as it trsvels 
54 km In one hour. Over a one degree square at 50° N about 400 cubic kilometers 
of air pass upward through the 200 meter level each hour in the area of strong 
convergence. 
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Figure 3 shows respectively the total mass of water vapor over each square 
meter of the sea surface in kg/m^ , the total liquid water in kg/m^ and the 
rainfall rate m mm/hr. The fronts In these figures are different from those 
in Fig. 2. T he analyses involved were independent and need to be reconciled 
on the basis of more data and consistency checks. 

These figures for REV 1212 provide valuable quantitive data on the properties 
of an extratropica L cyclone. The winds are shown to be transporting warm moist 
air by means of a concentrated southerly flow in advance of the cold front. 

The moisture is concentrated by this flow in a narrow band along the surface 
front. The band of high water vapor content undoubtedly would continue beyond 
the swath following along the frontal system. The highest concentrations of 
liquid water in Fig. 3(b) are along the front and in the area of convergence of 
Fig. 2(b). The precipitation is patchy in small cells with values over 4 iran/hr. 
Were the pattern in Fig. 3(c) to be shifted north north westward, without change 
in form, past a point at about 49° N, 130° W, at, say, 15 or 20 m/s and the 
rainfall rate integrated with time, the result would be an integrated iainfall 
^aph very similar to the one in Fig. 1(b). The cold front passage as in Fig. 1 
would also appear as a time history at this and many other points properly lo- 
cated in time. 


WIND STRESS AND WIND STRESS CURL 

There are many different equations that have been proposed to relate the 
wind at 10 m to the wind stress defined either by 

t/p - uj = - U V> (13) 

or by C 1Q * u2/ll(10)“ (14) 

as, i ' r example, Large and Pond (1981). A completely new one based on data 
from ttu^e different sources was used in this study. Lots of questions in this 
problem area remain unsolved, but if the superobservation represents the synop- 
tic scale wine at 19.5 meters, it can be related to the wind at 10 m by (15). 

0(19.5) = sJ ( 10) 4- (u*(U(10))£n(l9,5)/10))/K (15) 

and by, in meters per second, 

u* = 10- 3 (2.717(U(10)) + 0.142 (U(10) ) 2 + 0.0761(U(10)) 3 ) (16) 

For selected values of U (1 9 . 5 ) Table 1 can be obtained as examples and as de- 
tailed a table as desired can be constructed. There is a one to one relation- 
ship between the isotach values in Fig. 2(b) and u| values from the following 
table and the above equations except that u2 is strongly concentrated in areas 
of high winds. As the wind varies by a factor of 5 the stress varies by a 
factor of 45. 

2 2 

TABLE 1 Values of U(19.5), U(10), u^ and u^ in Meters per Second Except u^. 


U(19.5) 

U(10) 

2 

u * 

u * 

5 

4.7 

0.024 

0.156 

10 

9.5 

0.103 

0.321 

15 

14.4 

0.281 

0,531 

20 

18.7 

0.601 

0.775 

25 

23.2 

1.100 

1.049 


3 *+* 

For an air density of 1.25 kg/m , the vector value of x can be computed along 
with corresponding parallel and normal estimates of sampling variability. 

The finite difference estimate of tne curl is given by (17) in units of 
((Nevrtons per square meter) per meter). 
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Fig. 4 shows the field of the curl of the 
wind stress within the swath for REV 1212. 

For the four REVS of Pierson, et al. (1983), 
this figure is typical. Both high wind stress 
values and high values of the curl are con- 
centrated over relatively small areas of the 
ocean . 

At the synoptic scale, the wind stress is 
large over relatively small areas and small 
over relatively large areas. Areas of large 
wind stress curl are concentrated over small 
areas. The erratic day to day motions of 
drogued buoys as studied by Kirwan, et al. 
(1978) in the North Pacific may eventually be 
explainable by such data. 


DISCUSSION 

The time is not too f av off relative to a 
time scale beginning around 1900 and acceler- 
ating since SI ylab in 1973 and SEASAT in 1978 
to plans for several sc atterometers, radi- 
ometers and altimeters in the late 1980 , s when 
analyses such as this one over the oceans will 
be global and used to prepare routine weather 
forecasts and to understand the ocean circu- 
lation. There are problems of four dimen- 
sional space time data assimilation still to 
be solved plus the need for better atmospher- 
ic models and improved representations for 
turbulence. It will not be possible to beg 
the question by arguing that the initial value 
specification is not good enough to yield a 
correct forecast. It is already clear from the analysis of SEASAT products that 
the specification of wind fields over the oceans can, at times, be grossly in 
error when only conventional data are used. 



Fig. 4 The Curl of the Wind 
Stress for REV 1212. Values 
are in Newtons per Square Meter 
per Meter. 
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